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1 Introduction 

Let us recall the classical question by Fermi: "What plays the part of a hydrogen atom in this 
problem?" As to nuclear physics, certainly, its hydrogen atom is the deuteron. 

It is surprising how many nontrivial problems related to the deuteron can be solved by 
means of sufficiently simple, sometimes truly back-of-the-envelope analytical calculations. The 
reason is essentially that the deuteron binding energy, e = 2.23 MeV, is anomalously small on 
the nuclear scale. As a result of it, the deuteron wave function decreases very slowly beyond the 
range of nuclear forces. In the exponential asymptotics e~'^^ at long distances the parameter k 
is small: k = JmpE = 45.7 MeV {rrip is the proton mass). Correspondingly, the typical length 
at which the wave function falls down, is large: k^^ = 4.3 fm, it is much larger than the range 
of nuclear forces tq ~ 1 fm. This allows one to use for the deuteron wave function the so-called 
zero-range (of nuclear forces) approximation where this wave function is approximated by its 
asymptotic expression 




V'. = J^— m 



Due to the coefficient y^fi;/27r, this expression is normalized correctly: J dvtp'^ = 1. 

In the present review some results are considered which were obtained relatively recently 
basing on the zero-range approximation. Being related to the deuteron, these problems belong 
nevertheless to somewhat different fields: to the traditional nuclear physics, and to the problem 
of nonconservation of parity and time-reversal invariance. 

It is convenient however to discuss at first briefly few classical problems belonging already 
to textbooks or simply to physical folklore. The deuteron ground state contains in line with 
^5*1 a small admixture of ^Di, which can be in many cases neglected. The wave function of the 
triplet ^Si state of the continuous spectrum of the system neutron-proton can be written in the 
low-energy limit as 

^st = l- -, (2) 

r 

where at = 5.42 fm is the triplet scattering length. The analogous expression for the singlet 

^5*0 function of the continuous spectrum is 

i^ss = l- ^. (3) 

r 

The singlet scattering length is negative and very large: a^ = —23.7 fm. The subscript S at 
the wave function denotes here and below an S-wave, the subscripts t and s denote triplet 
and singlet states, respectively. Somewhat more precise expressions for the wave functions of 
5'-states of continuous spectrum, taking into account the relative momentum p of scattering 
particles (it is small as compared to l/rg, but not as compared to 1/at.s), can be conveniently 
written as |Q 

sin pr at,s e^^'' .^. 

i^st,Ss = -r-—- • (4) 

pr 1 + tpat^s r 

The simplest of the mentioned classical results is the relation between the deuteron binding 
energy and the triplet scattering length. From the orthogonality of the triplet wave functions 
(^ and (1^) which correspond to different energies, one obtains immediately 

at = -, or at = ^ (5) 

K y/rnpE 

It is valid within ~ 20%. 



Let us consider now the cross-section of the deuteron photodisintegration '-jd -^ np. It is 
natural to start with the contribution of El transition B R|. The dipole matrix element is 



'^l^l'^)=-4\/T^^rTx^- (6) 

V 1 — KTt (k^ + p^)^ 
The corresponding contribution to the cross-section is 

^ Svra Kj/ 

'^^^ 3 (l-fi:ri)(fi:2+p2)3- y') 

The result of a straightforward calculation with wave function (|I]) is augmented here with the 
correction factor (1 — /tr^)^^/^ in matrix element (^, and correspondingly, by (1 — nvt)^^ in 
cross-section (^; here rj = 1.76 fm is the so-called effective interaction radius. The correction 
is not so small, 1 — nrt = 0.59. Its origin is as follows. Matrix element (^) is dominated by large 
distances, and at these distances the account for the effective radius modifies the asymptotics of 
the deuteron wave function just in this way 0, |[. As usual, the cross-section of a El transition 
is proportional at the threshold to p^. 

Meanwhile, the contribution of Ml transition into the deuteron photodisintegration cross- 
section decreases in the threshold region as p^ only, and therefore dominates very close to the 
threshold |^, |^. Let us consider this contribution. Since the orbital angular momentum of 
the nucleons in the deuteron equals to zero (we neglect here the admixture of D-wave), the 
magnetic moment operator reduces here to a purely spin one: 

e 
M = [fiptTp + /i„cr„). 

Here a-p and cr„ are the spin operators of the proton and neutron, respectively, and fip = 2.79 
and fin = —1.91 are their magnetic moments. Due to the same orthogonality of the wave 
functions of the deuteron and '^5'i state of positive energy, the Ml transition goes into the 
singlet ^5*0 state of the continuous spectrum. A simple calculation |^, Q results in 






_ ^/iLt i%p yfj^p — ^n) \^ — t^^s) /qv. 



and correspondingly, 

2'na Kp (/ip — yUn)^(l — I^Oi, 

Though this contribution to the cross-section is enhanced by large numerical factors, fip — Hn = 
4.7 and 1 — nas = 6.5, it dominates only within 0.2 MeV above the threshold. 

After this, rather lengthy. Introduction we go over to the main part of the review. 

2 The n;?- scattering amplitude 

and the deuteron quadrupole moment 

Just as the triplet scattering length at is related to the asymptotic behaviour of the deuteron 
^5*1 wave function at r — >■ cxd, i.e., to the parameter k (see (^)), the limiting threshold values of 
the spin-dependent invariant amplitudes is related to the asymptotic value of the ^Di admixture 
in the deuteron wave function |^. 



We start with the standard expression for the np scattering amphtude, vahd under the 
assumptions of P and T invariance, as well as of the charge independence of nuclear forces 

(see 0): 

/(p'>P) = a + &(crin)(cr2n) + c(cri + cr2)n+(5f + /i)(crin+)(cr2n+) + (5f-/i)(crin_)(cr2n_). (10) 

Here p p' are the nucleon relative momenta in the center-of-mass system, initial and final, 
respectively; cri and (T2 are their spin operators; 

p ± p' p X p' 

n± = 1 — \ — TT ) n 



|P i P'l |P ^ P'l 

We get rid of the structure (crin)(cr2n) in ( PUJ ) by means of the obvious relationship 
(crin)(cr2n) = (Ti(T2 — (crin_(_)(cr2n+) — (crin_)(cr2n_). Then, we use the momenta p and 
p' instead of the unit vectors n±. And at last, we go over from the operators cri and (T2 to 
the total spin operator S = l/2(cri + 0-2) since we are interested in the triplet scattering. As a 
result, the triplet np scattering amplitude becomes 

/(p', p) = -a, + i^ |c'S[p' xp]+g, (S(p' + p))' + 92 (S(p' - p))^ , (11) 

Here the triplet scattering length is related to the parameters of formula (|l^) as follows: at = 
— (a + 6); the constants 

, 2cmp {g -h + h)ml (g -b- h)m^ 

^ ^ p^sinO ' ^' " 2p2cos2(^/2) ' ^' " 2p^sm\e/2) ' 

as well as at, are independent of the scattering angle 6. Expression (|Il|) is in fact an expansion 
of the scattering amplitude in momenta, where one retains in line with at only those higher 
order terms which depend on the spin 5*. 

We will go over now to the construction of the deuteron wave function at large r, which 
includes the D wave admixture. To this end, we find at first the effective ^-function poten- 
tial (pseudopotential) which reproduces in the Born approximation the scattering amplitude 
(ini) (see [H). With the correspondence 

p'iP'j -^ -ViV,(5(r) ,p[p, -^ -Vi6{T)V,,PiPj -^ -6{r)ViVj , 

we obtain 

Att c' 1 

f/(r) = — {at 5{r) + — e^,kS^ V, (5(r) V^ + -^% [{g, + (72) (V.V,<5(r) + 5(r) V.V,) (12) 
nip nip 2m^ 

+ (^i-^2)(V,5(r)V, + V/(r)V,)]}, 
where 



^ij — ^i^j 1 Ojbi dij. 



4 
3 

In the wave function of the scattering problem the outgoing spherical wave can be presented 
as 



m 



p 



gjp|r-r'| 



4^ / dr'j:—:jjU{r')Mr'). (13) 



We assume that the nonperturbed solution \l/o is an S wave one. Then expression ([13| ) reduces 
to 



By the analytical continuation of wave function ( |14|) to the point p = in, which corresponds to 



the bound state, we arrive at the following expression for the deuteron wave function at large 
distances: ^ 



The expectation value of the quadrupole moment tensor in this state is 

4 atm^ 

Now the deuteron quadrupole moment is 

Q = Q..\s.=i= 1^. (17) 

Let us note that, since the principal contribution to Q is given by r ~ 1/k, it does not make 
sense to introduce the correction factors (1 — Krt)~^''^ and (1 — nrt)''^ into expressions (|I5|) 
(p!6D, respectively. 

The analysis [P, P] of the experimental data on the elastic dp scattering and on the stripping 
reaction results in 

^i + (72 = 67fm. (18) 

The corresponding value of the deuteron quadrupole moment is 

Qi = 0.18 fm2. (19) 

Somewhat larger values follow from the phase analysis |T^ of the np scattering: 

Qi = (0.20 - 0.24) fm^. (20) 



At least the last value, (pO|), is not so far from the result of the direct experimental measurement 
of the deuteron quadrupole moment, 

Q = 0.286 fm^. (21) 



Of course, the accuracy of the direct experimental result (|2lD for the deuteron quadrupole 
moment is much higher then the accuracy both of the approximations made for the derivation 
of formula ([T7|) and of the phase analysis of the mentioned elastic processes and the stripping 
reaction. Therefore, the relations obtained, besides their theoretical interest, may be useful 
mainly for checks of phenomenological descriptions of these reactions. 

3 Circular polarization of 7-quanta 

in the reaction np -^ d^ with polarized neutrons 

Just as the photodisintegration of the deuteron at the threshold goes into the ^Sq wave of 
continuous spectrum, the radiative capture of thermal neutrons np — > d'j proceeds from the 



same ^Sq state via an Ml transition. Naturally, in this approximation, when the initial ^Sq state 
is completely spherically-symmetric, the polarization of an initial particle cannot be transferred 
to a final particle. However, due to the D wave admixture in the deuteron wave function and in 
the incoming wave, the Ml transition proceeds also from the triplet initial state. Besides, the 
corrections to the Ml operator responsible for the nonadditivity of nucleon magnetic moments 
in the deuteron also allow for a magnetic dipole transition from the triplet initial state. And at 
last, due to the D wave admixture an E2 transition becomes possible as well. All these effects 
are fairly small, but their investigation gives us an information on some subtle details of the 
np interactions at low energies. 

The circular polarization of the 7-quanta emitted in the radiative capture of polarized 
thermal neutrons by unpolarized protons np -^ d'~f was for the first time measured in [jTT| . This 
result, 

P^ = -(2.90 ±0.87) X 10"^ 

was improved essentially in the next work of the same group [0 : 

P-y = -(1.5 ±0.3) X 10"^ 



Theoretically the problem has been considered in articles [13-1^ 
essentially follow works [|^, |l^ . 



(22) 
In the present review we 



We start with the magnetic moment operator 

e 

— — [Ip + fJ^pCTp + fXr. 



M 



The proton orbital angular momentum Ip is obviously related to the total orbital angular 
momentum L: Ip = 1/2L. The linear combination of L and crp„ can be presented as follows: 



- L ± fipCTp + /i„cr„ 



L±-o-p 



1 



1 



+ 2(/^p 



yUn)(crp - cr„) + - { fip + fin - 



[CTp + CTn] 



(23) 



_ 1 

The first term in the second line is nothing but one half of the total angular momentum I which, 
being an integral of motion, cannot cause any transition. The second term is responsible for 
the common Ml transition from the ^5*0 state of the continuous spectrum, i.e., for the process 
inverse to the deuteron photodisintegration at the threshold. The effective operator of this, 
principal Ml transition, or the coordinate matrix element of the second term, is 

fIn)v27TK, 



Mn 



2m„ 2 ^^^ 



K^ 



(24) 



As to the last term in the second line of (p3D, its matrix element is distinct from zero first 
of all due to the ^Di admixture to the wave functions of the deuteron and the ^Si state of the 
continuous spectrum, ipd and ipst, respectively. For calculating this matrix element we use the 
orthogonality condition for the complete radial wave functions of the both states: 



drr^RdoRto + Rd2Rt2) = 0. 



(25) 



Here we mark by the second subscripts, and 2, the radial functions of the S and D components 
of the initial and final states. With the account for (|25|) the discussed matrix element is found 
easily: 



lip + fin- 7T ) {f^p + (^n) 



2,/^P + /^n 2 



I 



drr^R^oR 



d2J^t2 



0. 



(26) 



The last radial integral is dominated by distances smaller than the range of nuclear forces. It is 



natural to assume |]T3| that at these distances Rd2 and Rt2 differ in normalization factors only, 
i.e., that 



R 



42 



-at 




2tx 

Rd2- 

K 



Then matrix element 



reduces to 
3 



where Pd is the D wave weight in the deuteron. 

Let us recall now the known expression for the deuteron magnetic moment ^d- 

l^d = l^p + l^n - ^ifJ'p + fJ'n- ^jPd + '^IJ' 



(27) 



(2^ 



(numerically, Hd = 0.8574). In line with the correction for the D wave admixture (the term 
with Pd), it includes the contribution of relativistic corrections A/x. Naturally, the matrix 
elements of corresponding operators are also dominated by short distances. Let us assume 
again that at these distances the wave function of the ^5*1 state of continuous spectrum differs 
from the deuteron one in the overall factor only. Then the coordinate matrix element of the 
corresponding operator of the triplet-triplet transition is 



-at 




K 



Arising in this way total effective operator of the Ml triplet-triplet transition is 



M 



2rar 



or, in virtue of (p! 



M 



2 f/^P + /^n- 2 1 ^'i 



(yUp + /in - ^^d)at 



Afi 



2tt 



2m, 




at 



I. 




I, 



K 



K 



(29) 



(30) 



A standard calculation, taking into account expression ( p^ ) for the principal transition, 
leads to the following result for the Ml contribution to the degree of circular polarization: 



P^(Ml) = - 



fi'P + fJ'n- fJ'd i^at 



Up fi,^ 



1 



-0.92 X 10" 



KCto 



(31) 



One more contribution to the degree of circular polarization of photons is due to the 
quadrupole transition. The E2 operator equals here 

O 

(let us recall the relation r^ = r/2 between the proton coordinate r^ and the argument r of 
the wave function). Quite standard calculation with wave functions (0) and (plSf ) leads to the 
following result for the quadrupole contribution to the circular polarization: 



P,{E2) 



K-^{9l+92) 



15 ml (1 - Ka^) {jJ^p - jJ^n) 



2 atK^Q 

5(1- Kas) i^Ip - fin) 



-0.25 X 10" 



(32) 



Our final result 

P., = P^(Ml) + P^{E2) = -1.2 X 10-3 (33) 

and the experimental result (^) are in a reasonable agreement. There is also a good agreement 
between (^) and the final results of recent calculations [|17|, |1^), though for P^(Ml) and 
P^{E2), taken separately, our and their results differ considerably. 

4 P- and T-odd electromagnetic moments of deuteron 

4.1 Generalities on anapole moments 

A notion of anapole moment (AM) was introduced by V.G. Vaks (who was then a graduate stu- 
dent of A.B. Migdal) and Ya.B. Zel'dovich |T^. AM is a special electromagnetic characteristic 



of a system where parity is not conserved. 

The peculiarity of AM is in particular that the interaction of a charged probe particle with 
an anapole moment is of a contact nature (for a more detailed discussion see, for instance, |]20(| ). 
Therefore, for example, the interaction of the electron with the nucleon AM, being on the order 
of aG {G is the Fermi weak interaction constant), cannot be distinguished in general case from 
other radiative corrections to the weak electron-nucleon interaction. And in a gauge theory of 
electroweak interactions only the total scattering amplitude, i.e., the sum of all diagrams on the 
order of aG, is gauge- invariant, independent of the gauge choice for the Green's functions of 
heavy vector bosons. Thus, generally speaking, the AM of an elementary particle or a nucleus 
is not gauge-invariant notion, and therefore has no direct physical meaning. However, there are 
special cases where the anapole moment has a real physical meaning. In heavy nuclei the AM is 
enhanced ~ A^^^ [^ {A is the atomic number), as distinct from common radiative corrections^. 
In fact, the nuclear AM of ^^^Cs was discovered and measured with good accuracy in an atomic 
experiment |2^. The result of this experiment is in a reasonable quantitative agreement with 



the theoretical predictions, starting with I^T], Q. 

There is one more object, the deuteron, whose anapole moment could make sense for a 
sufficiently large P odd vrNN constant [^. The problem of the deuteron AM was considered 
phenomenologically in [1,23-25] . Recently the deuteron AM induced by the P-odd pion exchange 
was calculated in pS) (see also p7|). 



4.2 Nucleon anapole moment 

It is convenient to start the discussion with the nucleon AM in the limit mT^ -^ 0. It was 

calculated in 1980 by A.I. Vainshtein and the author. The result is the same for the proton 

and neutron: 

egg f 6 m^ mp\ 
ap = a„ = - — 1 In ^^ <t; 34 

l/rripTTlTr V ^ "^p fT^TT J 

we assume e > 0. Being the only contribution to the nucleon AM, which is singular in rriT^, the 
result (|3^) is gauge-invariant. In this respect, it has a physical meaning. 

Unfortunately, despite the singularity in m.,r, the corresponding contribution to the electron- 
nucleon scattering amplitude is small numerically as compared to other radiative corrections 



^It means, by the way, that there is an intrinsic hmit for the relative accuracy, ~ A^^^, with which the AM 
of a heavy nucleus can be defined at all. In the case of ^^'^Cs, which is of the experimental interest, this limiting 
accuracy is about 4%. 



to the weak scattering amplitude. Indeed, the radiative corrections to the effective constants 
C2p,n of the proton and neutron axial neutral-current operators G/\^C2p,nCrp,n are [^ 

C^p = -0.032 ± 0.030 , C;^ = 0.018 ± 0.030 . (35) 

In the same units G/\/2, the effective axial constants induced by the electromagnetic interaction 
with the proton and neutron anapole moments (0), is 

C;_„ = -aaN{\e\G/V2)-^ = 0.07 x 10^^. 

At the "best value" g = 3.3 x 10^^ (strongly supported by the experimental result for the ^^^Cs 
anapole moment) we obtain 

Cp% = 0.002. (36) 

With this value being much less than both central points and error bars in (|35|), the notion of 
the nucleon AM practically has no physical meaning. This is why the result (|3^ ) was never 
published by the authors. It is quoted in book [|^ (without the logarithmic term) just as a 
theoretical curiosity. This result was obtained also in ||2^, the logarithmic term in the nucleon 
AM is discussed in |^ . 



Let us note, that as pointed out in ||311|, a P-odd ttttNN interaction also generates a contribu- 
tion ~ InrTijr to the nucleon AM. Only purely theoretical estimates are known for the constants 
of this P-odd vrvrNN interaction. According to the estimates considered by the authors of |^| 



as relatively reliable, the contribution of the P-odd vrvrNN interaction to the nucleon AM is 
about an order of magnitude smaller than (|34D . 

However, the situation with the deuteron AM is quite different. Not only the proton and 
neutron AMs add up here. The isoscalar part of the radiative corrections is much smaller than 
the individual contributions Cgp and CJ^, and is calculated with much better accuracy p8| : 

Cm = C2p + CL = 0.014 ± 0.0030 . (37) 

Moreover, there are new large contributions ~ rra"^ to the deuteron AM, induced by the P-odd 
TT-meson exchange. We go over now to their consideration. 

4.3 P-odd TT-meson exchange 

The Lagrangians of the strong vrNN interaction and of the weak P-odd one, Lg and L^, respec- 
tively, are well-known: 

Ls = g[V2 [pi'j^n 7r+ + ni'j^p vr") + ( pi'j^p - ni'j^n) 7r°] ; (38) 

L^j = g V2i{pn7r^ — npn~). (39) 

Our convention for 75 is 

7.= (J-J); (40) 

the relation between our P-odd ttNN constant g and the commonly used one h^^lfjsf is gV^ = 
h-lfAT- Our sign convention for the coupling constants is standard: g = 13.45, and (7 > for 



the range of values discussed in |^ 



The resulting effective nonrelativistic Hamiltonian of the P-odd nucleon-nucleon interaction 
in the deuteron due to the pion exchange is in the momentum representation 

^(q) = — -^^(Nln.N,) (nIt^^n^). (41) 



Here q = pi — Pi = — (P2 — P2) = Pn ~ Pp = ~(Pp ~ Pn)- This P-odd interaction conserves the 
total spin 

(and does not conserve the isotopic spin). Thus in our problem it mixes the ^Pi state of the 
continuous spectrum and the deuteron ground state ^Si. Let us note that the P-odd interaction 
( PD which interchanges the proton and neutron, when applied to the initial state aU ri)al{r 2) \0) 



transforms it into a|^(ri)a^(r2)|0) = — a]^(r2)a|j(ri)|0). On the other hand, the coordinate wave 
function of the admixed ^Pi state is proportional to the relative coordinate r, which we define as 
Tp — r„. Therefore, it also changes sign under the permutation p ^^ n. Thus, for the deuteron 
the P-odd potential can be written in the coordinate representation as a simple function of 
r = Fp — r„ without any indication of the isotopic variables: 

V{r) = -^l- {-^I . V) '.^^tUhlll . (42) 

The corresponding imaginary mixing matrix elements are related as follows: 



the overall sign in (^2]) corresponds to {^Pi\VfSi). 

The weak interaction (|42|) generates a contact current j'^. To obtain an explicit expression 
for it, we have to consider the P-odd interaction V in the presence of the electromagnetic 
field. Its including modifies the proton momentum: p — ^ p — eA, which results in the shift 
q ^ q + eA in interaction (pt). Then in the momentum representation the contact current is 



J (q) 



dV{q) d 2gg I(q + eA) 



dA dA nip ml + (q + eA)^ 

2egg ( I 2q(Iq) 



rup \ml + q^ {ml + q^y ' 

In the last expression we have neglected the dependence on A. In the coordinate representation 
the contact current is 

rir) = ^ r(IV) ^^ . (44) 

27Tm,p r 

4.4 Calculation of the deuteron anapole moment 

Let us discuss at first a general structure of the deuteron AM generated by a P-odd np interac- 
tion, which conserves the total spin, assuming only that the deuteron is a pure ^5*1 state bound 
by a spherically symmetric potential. We follow here essentially the arguments applied in p3 



see also book [EOl) to the problem of a single proton in a spherically symmetric potential. In 



this case of a single proton the AM operator is |23 



— Uprxa--[f,r]^+—rx[rxf]. (45) 



In the case of the deuteron this formula generalizes to 



ad = -^ < r X {fipo-p - ^„o-„) - - [P, r] ^ + - r X [r X j'^]. (46) 



Both operators (|45| ) and (^) are orthogonal to r (neither of them commutes with r, so the 



orthogonahty means here that ar + ra = 0). Therefore the contact current (PI), generated by 
the P-odd pion exchange and directed along r, does not contribute to the nuclear AM. 

Let us present the wave function of the deuteron ^Si state as ipQ{r)xt- If the weak interaction 
conserves the total spin I, the P-odd admixture ^Pi can be written as 5V'i(r) = i{Ir/r)i{ji{r)xt 
(both radial wave functions, ipo{^) and ^i(r), are spherically symmetric). Simple calculations 
demonstrate that the expectation value of operator (|^) in the state with wave function [iljQ{r) + 
i{lr/r)ilji{r)]xt is (in the absence of a contact current contribution) 

ad = - — inp- fin- -jl dYipQ{r)ripi{r) . (47) 

So, under the assumptions made, the deuteron AM should depend on the universal combination 

yUp - yUn - 1/3. 

The real calculations we perform in the zero-range approximation where ^^^{r) = ^d(r) 
(see (IID), and the wave functions of P states are free. Moreover, when using the stationary 
perturbation theory, we can choose plane waves as intermediate states we are summing over, 
since perturbation V^(r) (see (^21) ) will select by itself P states from plane waves. The correction 
to the wave function is 

/rlV p'kr f 

(2^373117^ /*'^-"^M'^o(0. (48) 

Rather lengthy calculation leads to the following expression for the matrix element of the 
radius- vector: 

rfrV'd(r)r5^i(r) = -z-^^-^^t^TT, (49) 

where ^ = n/mT^ = 0.32. With this matrix element and operator (|46| ) one obtains easily the 
following expression for the deuteron AM: 

in accordance with the general formula (|47|). 

In fact, the range 1/m.^ of the P-odd interaction (^) is quite comparable to the range of 
the usual nuclear forces. Therefore, it is rather dangerous to use here blindly the naive wave 
function (|l]). Still, numerical calculations with a model wave function, which has more realistic 
properties, lead to the result which differs from (^) no more than by 20%. As to other sources 
of P- violation, different from the pion exchange, it can be demonstrated that their contribution 
to the deuteron AM does not exceed 5% of ( ^Ol) (at the "best value" g). 

It looks reasonable to combine the potential contribution ( |50D with the sum of the proton 
and neutron anapole moments (|3^) . In this way we arrive at the final result for the deuteron 
AM: 

ad= --^^^[0.49(/ip-/i„-l/3) + 0.46] I = -2.60-^^^1. (51) 

This result includes all contributions to the P-odd amplitude of erf-scattering, which are singular 
in m^r, and thus is gauge- invariant, physically sensible. 

Now let us compare the contribution of the found AM (^) with the P-odd erf-scattering 
amplitude due to usual radiative corrections, which are nonsingular in m.,r- At the "best value" 
g = 3.3 X 10"^ and the estimate of 20% for the accuracy we obtain 

C^a = 0.014 ± 0.003. (52) 



This number is quite comparable with the contribution (|37|) of usual radiative correction. Taken 
together, they constitute 

C^d = Cl^ + C^d = 0-028 ± 0.005 . (53) 



One more contribution to C2d is due to the admixture of strange quarks in nucleons |]33[ . The 
magnitude of this contribution is both extremely interesting and highly uncertain. 

To measure experimentally the constant C2d is very difficult. However, with such a good 
accuracy of the theoretical prediction (p3D , this experiment would be a source of valuable 
information on the P-odd vrNN constant and on the strange quark admixture in nucleons. 

4.5 P-odd, T-odd electromagnetic moments of the deuteron 

P-odd, T-odd multipoles of the deuteron, electric dipole and magnetic quadrupole, were previ- 



ously considered phenomenologically in |34]. Then in |2^ these multipoles have been calculated 
in the same approach as the anapole moment. 

Three independent P-odd, T-odd effective vrNN Lagrangians are conveniently classified by 
their isotopic properties: 

AT = 0. Lq = gQ[V2{pmT^ + npTT~) +{pp-nn)7T% (54) 

|Ar| = 1. Li = giNNn^ = gi{pp + nn)7i°); (55) 

|AT|=2. L2 = g2iNTN7z - 3Nt^Ntt°) 

= g2 [V2(j)n7r^np7c^) — 2{pp — nn ) vr*') ]. (56) 

Since the possible values of the isotopic spin for two nucleons is T = 0, 1 only, the last interac- 
tion, with I AT I = 2, is not operative in our approach. 



The effective P-odd, T-odd np interaction is derived in the same way as (^21) . In the 
momentum representation it looks as follows: 

W{r) = -^—[{3go-gi)(Tp - {3go + gi) (Tn]V ^^ . (57) 

The calculation of the deuteron electric dipole moment d^^, i.e., of the er^ = er/2 expectation 
value, is performed analogously to that of the AM. The result is 

d = -^gg^ / + ^, T. (58) 

127rm^ (1 + 20^ ^ ' 

The magnetic quadrupole moment (MQM) operator is expressed through the current density 
j as follows (see, e.g., [EI 



^^^mn ym^nrs ' 'n^mrsj'rjs- \^'^) 

This expression transforms to 



M^n = — <J 3/i 



......„-?(.. 



+ 2q{rmln + Tnlm) \ ■ (60) 



Here fi is the total magnetic moment of the particle, q is its charge in the units of e. The 
magnetic quadrupole moment is the expectation value A4 of the operator M^z in the state with 
the maximum spin projection I^ = I. 



In our case, due to the spherical symmetry of the deuteron nonperturbed wave function, 
the orbital contribution to Mmn vanishes. The contact current here also is directed along r and 
therefore does not contribute to MQM. Thus the deuteron MQM arises due to the spin term 
in (ppp. It equals 

M = -I^^^(I^[(3^o + ^i)/..+ (3^o-^i)/.n)]. (61) 

Let us note in conclusion that the possibility to search for the deuteron dipole moment at a 
storage ring with polarized nuclei is discussed at present by experimentalists rather seriously. 

*** 
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